Holonomic quantum computation in the presence of decoherence 
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We present a scheme to study non-abelian adiabatic holonomies for open Markovian systems. As 
an application of our framework, we analyze the robustness of holonomic quantum computation 
against decoherence. We pinpoint the sources of error that must be corrected to achieve a geometric 
implementation of quantum computation completely resilient to Markovian decoherence. 
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Adiabatic holonomies in quantum mechanics are uni- 
tary transformations generated by slowly changing the 
Hamiltonian of a quantum system through a set of pa- 
rameters in a cyclic fashion. The parameters define a 
manifold and the holonomies depend only on the path 
followed by the system in this space. In case the Hamil- 
tonian is non-degenerate, the holononry is an abelian 
phase better known as Berry's phase If the Hamil- 
tonian is degenerate, the holonomy is a non-abelian gen- 
eralization of the phase which induces transitions among 
states belonging to the degenerate subspace 0- Geo- 
metric phases, i.e.. Berry's phase and its various gen- 
eralizations, were extensively studied in the 1980's 3] 
and recently became, once more, a fashionable subject 
since it was proposed that quantum computation (QC) 
could be implemented in a fault-tolerant way using geo- 
metric transformations The most general scheme 
of this kind is called holonomic quantum computation 
(HQC) and it involves constructing a universal set of 
quantum gates by making the system acquire a succes- 
sion of abelian and non-abelian geometric phases. Since 
holonomies depend only on geometric properties of the 
parameter space, HQC is robust, by construction, to any 
kind of error that leaves these properties invariant. This 
is considered to be the strongest advantage of the scheme. 

In any realistic implementation of quantum computa- 
tion one must consider that, in practice, quantum sys- 
tems are never isolated. The interaction with the envi- 
ronment effects that initially pure quantum states decay 
into statistical mixtures (called mixed states) through a 
process known as decoherence It is known that this 
process is a severe limiting factor for quantum compu- 
tation. Therefore, in a pragmatic realization of HQC 
schemes, it is crucial to define holonomies for mixed 
states, and understand how they are generated under the 
presence of decoherence. So far, non-abelian holonomies 
have been only investigated for pure quantum states 
evolving under unitary transformations. It is now neces- 
sary to investigate to what extent the geometric proper- 
ties on which holonomies depend, are left invariant under 
quantum noise. 

Although claims that quantum gates of geometric ori- 
gin offer increased fault-tolerance can be found through- 
out the literature, a comprehensive analysis on the ro- 



bustness to errors of these schemes has not yet been pre- 
sented. Partial results, however, have been obtained: For 
both, abelian and non-abelian phases, some errors of clas- 
sical origin have been investigated [7| and the effect of 
quantum noise has been analyzed only for the abelian 
case in 8, 9, lOj. In particular, Ref. ^ uses the quan- 
tum jump approach to study the effects of decoherence in 
the adiabatic and non-adiabatic case. There the spin-1/2 
Berry phase was found to be robust to the depolarizing 
channel in the no-jump trajectory. Moreover, resilience 
to the phase flip error is guaranteed for any number of 
jumps. In case the particle is driven by a quantized field, 
Berry's phase is found to be robust to field decoherence 
to second order 10]. 

The robustness of the abelian phase to a number of de- 
coherence effects suggests that non-abelian phases might 
behave similarly. Motivated by this, in this Letter, 
we investigate the effects of decoherence on non-abelian 
holonomies under the Markovian assumption, where en- 
vironmental memory effects are negligible. We demon- 
strate the application of our general framework to the 
study of decoherence in HQC. Independently of the phys- 
ical implementation considered, we estimate the errors, 
produced by generating a universal set of holonomic gates 
in the presence of an environment. We pinpoint to which 
specific types of errors the scheme is fallible. These re- 
sults constitute the basis of a framework to perform uni- 
versal HQC resilient to Markovian decoherence. Deco- 
herence can produce two kinds of errors: those which 
take the system out of the degenerate subspace and those 
which take place within the subspace. The first kind of 
error can be eliminated by working in the ground state 
and considering a system where the energy gap with the 
first excited state is very large. Our analysis concerns 
the second kind of errors. We consider only adiabatic 
holonomies because the non-adiabatic case is not well 
understood. It is not clear how to separate, in gen- 
eral, dynamic from geometric evolution in non-adiabatic 
holonomies. 

Our work focuses on the applications of holonomies in 
quantum computation, but geometric phases are a fun- 
damental aspect of quantum mechanics with applications 
in many fields. They have manifestations that range 
from low to high energy physics, appearing in optical and 
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solid state systems, in molecular and atomic physics, and 
are at the heart of phenomena such as anomalies in gauge 
field theories, fractional statistics, Aharonov-Bohm and 
quantum Hall effect. The discussion presented here is 
relevant to any physical situation involving holonomics 
where decoherence plays a role. 

To introduce the concept of holonomy we consider a 
Hamiltonian H{Xq) = Hq describing an energy degener- 
ate N-dimensional quantum system. The Hilbert space 
has k subspaces with corresponding energy Ek ■ Our anal- 
ysis will be restricted to one of the degenerate subspaces 
with rt-fold degeneracy and energy E. The system is 
initially prepared in a state belonging to this subspace 
= IV'i(O)) and the Hamiltonian of the system varied 
through a set of parameters A on a control manifold A4 
in an adiabatic way. The adiabatic theorem states that 
if the Hamiltonian is varied sufficiently slowly with re- 
spect to any time-scale associated with the dynamics, the 
system will remain in a state belonging to the subspace 
corresponding to energy E. In other words, there is no 
level-crossing. The degeneracy structure of the Hamilto- 
nian is preserved throughout the evolution which we can 
write as H{X) = V{X)HoV^X), where V{X) is unitary If 
the Hamiltonian returns to its initial value after a time T, 
H{T) = Hq, describing a closed curve C in A^, the state 
Itp^) is mapped to e~*^-^t7c(A)|'!/'°). The transformation 
Uc{X) is called the holonomy and in the following, we 
will derive it. 

By changing the Hamiltonian through the set of pa- 
rameters A, the state of the system is parallel trans- 
ported in parameter space. A rule for parallel transport 
in the manifold, i.e. a connection, must be specified since 
there is no unique way of parallel transporting a vec- 
tor. A connection is provided by requiring that the state 
\ipaiX)) remains normalized through parallel transport, 
(^laf^lV-a) =0. 

In terms of the local reference basis of the degenerate 
subspace {|^a(A))}, the state at any point of the adia- 
batic path in parameter space is expressed as, |-0c((A)) = 
Uafs{X)\(j)f3{X)), where U is unitary and |V'q(A)) corre- 
sponds to the solution of the Schrodinger equation for 
the initial condition = |(/)q(0)). The parallel trans- 
port condition then reads, 

Uls^Pp', + ulsUa(^A^, = 0, (1) 

with Pi3j = {4>j\4>/3) and A^-y = (<?i'7l gf^l*/!/?), in the 
Wilczek-Zee notations. The matrix P is hermitian 
and in the case when the states of the local basis are 
orthogonal it is equal to the identity. In the presence of 
decoherence, the latter is not true in general. We ignored 
the constant term EP/^j which produces the global dy- 
namical phase e~*^"^. Solving for U when P is invertible, 
U^^U — —AP^^, and integrating over the closed path C 



we obtain the holonomy, 

Uc{X)=Ve-^cAP-\ (2) 

where P is the path ordering operator. In terms of 
the eigenstates of the initial Hamiltonian TJq, we have 
Ap^ = (V'7l^^^l''/'j3)- The dimensionality n of the holon- 
omy equals the degree of degeneracy of the eigenspace. 
Berry phase is the special case when the eigenspace is 
non-degenerate, and the unitary transformation is then 
one dimensional, i.e. a complex number. The holon- 
omy (0) depends only on the path followed in parameter 
space, and transforms, under gauge transformation g, as 
Uc{A) gUc{A)g^^ We point out that all previous 
considerations can be extended to the case of an open 
curve without any modifications since gauge invariance 
is not relevant in our context. 

To investigate the effects of decoherence in the non- 
abelian geometric evolution of states described above, we 
employ the quantum jump approach. The master equa- 
tion (h^i) 

1 1 " 

p=-[H{X),p]~-J2{LlLkP + pLlLk -2LkpLl} (3) 
* fe=i 

dictates, in the Markovian approximation, the evolution 
of the system described by the density operator p. The 
commutator generates, through the Hamiltonian H(X), 
the coherent part of the evolution and the second part 
represents the effect of the environment on the dynam- 
ics of system. The operators L^, are called Lindbladian, 
and by prescribing them, one models different decoher- 
ing processes. Equation Q is in general very difficult to 
solve but the quantum jump approach provides an inge- 
nious solution to this problem. Consider that the time 
evolution of the density matrix, for small time intervals 
At, can be written as 

n 

p(< + Ai)«^W^,p(i)W^t^ (4) 

k=0 

where the operators Wk are complete positive maps ful- 
filling the completeness relation J2^=o^k^k = 1- By 
setting Wo = 1 - iHAt and Wk = VAtLk with H a 
non-Hermitian effective Hamiltonian, 

n 

H^H~'-Y,Llu, (5) 

fe=i 

the dynamics of the system is approximated by dividing 
the total evolution time T into a sequence of discrete in- 
tervals At = j^. Wq and Wk are called the "no-jump" 
and jump operators respectively. According to Eq. Q), 
the state of the system, after any time step t„i — mAf, 
evolves into p{tm+i) — Wkp{tm)Wl (up to first order in 
At), with probability pk = TrWkp{tm)Wl. The time 
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evolution of the system is then calculated for a set of 
possible trajectories containing, each one of them, differ- 
ent numbers of jumps occurring at different times. Each 
trajectory is defined as a chain of states obtained by the 
action of a sequence of operators Wk on the initial state. 
For example, for an initial pure state {tpo), the (non- 
normalized) state of the system, after the m-th step along 
the i-th trajectory, is given by: 

m 

\i^^^)^l[^^{l)Ho), (6) 
1=1 

where stands for the l-th element of a sequence 
of indexes with values from . . .n. Each trajectory 
is represented by a discrete sequence of pure states 

{■00 7 ''l'o \ ■ ■ ■ 1 ^'n}- dynamics given by the master 

equation is recovered by summing incoherently all the 
states associated to each trajectory, and taking the con- 
tinuous hmit At — > 0. 

The no-jump trajectory corresponds to the case in 
which no decay occurs. The evolution of a quantum state 
along this trajectory is obtained by the repeated action of 
the operator Wq and taking the continuous limit N ^ oo. 
This yields a dynamics governed by the complex effective 
Hamiltonian H: 

^l\i,it))^H\m), \m)^\^o). (7) 

Since the Hamiltonian is non-hermitian, the correspond- 
ing eigenstates are non-orthogonal. 

We are now ready to consider the case in which a 
non-abelian holonomy is generated by the transforma- 
tion H{X) — V{X)HoV^{X) when the system's dynamics 
governed by © . Note that it is only possible to generate 
the holonomy when the interaction of the system with 
the environment is such that the degeneracy structure 
of the Hamiltonian is preserved. This is the case when 
K = X]fc=i ^jc-^k = ctl OT K = aH. The holonomy in the 
no-jump trajectory is then the same as the one acquired 
by an isolated system evolving under the same Hamilto- 
nian H. This is because the eigenstates of H coincide 
with the eigenstates of H, so that P is proportional to 
the identity and the connection in the no-jump trajec- 
tory (AP"i)(") is equal to A = V^V. The interaction 
with the environment only produces an overall visibility 
factor e°*/2 foj. k = al and e"^*/^ foj. ^ ^ ^fj xhe 

factor is small for low decoherence rates and in the sec- 
ond case, it is eliminated by working in the ground state. 
In other words, for sources of decoherence that satisfy 
these conditions, the holonomy is robust in the no-jump 
trajectory. 

Let us consider a more interesting case, the trajectory 
in which there is one jump Wi at Ai. During the adiabatic 
evolution, the system evolves under a no-jump trajectory 
from Ao to Ai, then the jump occurs instantaneously 
and the system continues to evolve by the transforma- 
tion V{X) until time T, corresponding to A/. Using the 



composition rule for holonomies, Ui = (Pe ^^)(Pe 
where 

To = r\AP-i)(") = A, (8) 
Jo Jo 

Ti = {AP-^f \ (9) 

Since the one-jump connection for W}Wi = ail is given 
by, 

{AP-'i^^ = mv\X)V{X)W,)Pi^\ (10) 
10;) = l^.er°|<^°), (11) 
Ppi = (0710)3) = 0*^/37, (12) 

the holonomy after the jump is 

Pe-ri=Pe-^-f'/^^^^-. (13) 

So in the one-jump trajectory Ui ~ Pe~°i"'^''"°^'Pe~'"°. 
This result can then be generalized to a trajectory for 
which n jumps occur 

n 

K = HP'^^"*'^'^"^^ '<"Pe-^" (14) 
1=1 

rf = [''"a. (15) 

Note that we can write, 
p,-^w:,>r:rw.,, ^ l^t^^ (^Pe-/.7-4) W^^.y (16) 

We now consider the application of our results in QC, 
where the idea of using quantum systems to perform 
computations more efficiently than classical computers 
is investigated ■ Geometric phases play an important 
role in this field since quantum gates can be implemented 
in a geometric way. Jones et. al. 3 proposed a scheme 
for generating quantum single and two-qubit phase gates 
using Berry's phase. 

The generalized scheme to implement QC by geometric 
means is HQC |5|. In this scheme the input information 
is encoded in the n base states of a given degenerate sub- 
space at H{Xo). By renormalizing one can choose the 
energy of the subspace E = 0. The gates are gener- 
ated by the holonomic evolutions described above (this 
includes abelian and non-abelian phases) which are pro- 
duced by slowly changing the initial Hamiltonian through 
the set of parameters A. After A completes a loop C in 
rooted at Aq, the initially prepared state lipm) in which 
the information was encoded is mapped to an output 
state \ipout) = UciX)\'ipin) where Uc{X) is the quantum 
gate. To perform a given gate it is necessary to find 
the path or succession of paths in parameter space which 
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yields such a gate. A quantum algorithm is built from a 
sequence of gates such that the final state corresponds to 
the solution of the computational problem. Any unitary 
transformation, thus any algorithm, can be approximated 
by a strategical succession of closed paths in the control 
manifold M.. As long as the adiabatic condition holds, 
the computation does not depend on the rate at which 
the control loops are driven. The gates performed in such 
a way depend only on geometric properties of the param- 
eter space. Thus, any computational error which is path- 
preserving does not change the transformation. Hence, 
errors such as fluctuations in the driving parameters and 
systematic errors are automatically avoided. However, 
the environment might induce some fluctuations which 
are not path-preserving; it is then important to quantify 
them. 

Investigating the effects of decoherence in HQC be- 
comes crucial for the scheme to be realizable. Our quan- 
tum jump formalism can be used to investigate decoher- 
ence in a universal set of gates generated by holonomies. 
By universal we mean that any unitary evolution can 
be constructed from this set of elementary gates [llj . 
The set consists of two 1-qubit gates of the form Ui = 
giSio-i .yj,]^g]-g a^j-g Pauli matrices, and one 2-qubit gate 
[/g ^ (,i<l><y^®<y, _ For example, Ui = e^'^^''\ U2 = e**^'^^ 
and C/3 = e^'P'^^'^'^^ are universal since any transforma- 
tion in S'C/(4) can be generated from them. Consider 
that these gates are performed by holonomic evolution 
as described in but the system interacts with the 
environment. Errors within the subspace for one qubit 
gates are proportional to cr_ , cr+ and (7^, with i = 1, 2, 3. 
For Li = \f(xai and a the decoherence rate of the error, 
we obtain L\Li = al. Thus the universal set of gates is 
robust to these errors in the no-jump trajectory. When 
there is a jump in the trajectory, it is easy to see from 
eq. H16() that the gate Ui is robust to errors Lj when [4 



The analysis of the two-qubit gates follows directly 
from our previous conclusions. As an example we con- 
sider errors of the form Oi ® Uj. For simplicity let us 



consider the gate 



^ any other gate can be 



i — j and for errors with j ^ i one obtains a change in 
sign (since crjtTiCrj = —<Ji). For example, if one jump at 
occurs at Ai, the gate Ui becomes Ui — expi(6'i — 02)o'i 
where 

91= 02= n. (17) 

Here corresponds to the solid angle subtended by the 
evolution of the system in parameter space. The angle 
9 = 61 + 02 corresponds to the gate without errors. Gen- 
eralizing to n errors of the type Lj we find that Ui(n) = 

gie,^. with 9e = Em(-l)'"6'm ^ud 9,n = Q. When 

errors a± occur, the degeneracy structure of the Hamilto- 
nian is preserved only for initial Hamiltonians which are 
proportional to 0-3 since (J±azfz = i(l ± 0-3). In this case 
the holonomy is robust in the no-jump trajectory. For a 
single jump cr±, P is not invertible but we can calculate 
the gates from eq. H Unfortunately, the gates are com- 
pletely lost, e*^'^' becomes a U{1) phase after one jump. 
This is the main source of error that must be corrected. 



analyzed in the same way. The gate is robust to any 
number of jumps of the type (Ti (g) cti, cti (g) 1, 1 (8) ai and 
(T2 (8) (T2 but changes in signs occur for 1 (g) 0-2, fT2 1, 
(Ti (g) 0-2 and 0-2 (gi cTi . Nevertheless in general the most 
common non local error is ai^ ai, to which the 2-qubits 
gate is robust. 

We have presented a general scheme to study the ef- 
fects of Markovian decoherence in the generation of non- 
abelian adiabatic holonomies. We applied it to analyze 
the effects of the environment on a universal set of holo- 
nomic quantum gates. A scheme for quantum computa- 
tion completely robust to Markovian decoherence can be 
constructed by using HQC and finding a way to correct 
for the errors pointed out in this Letter by other means. 
This could possibly be achieved using error correction 
techniques. An example of the use of geometric phases 
to produce an error correcting code is explored in [l^ . 
We are currently working on the details of such a scheme. 
Since holonomies are a main ingredient in loop quantum 
gravity it is likely our scheme can be applied to study 
decoherence in this field. 
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